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Abstract 

We consider the possibility that there exist sterile neutrinos which are closely degenerate 
in mass with the active neutrinos and mixed with them. We investigate the effects of this 
kind of active-sterile neutrino mixing on the composition of supernova neutrino flux at the 
Earth. If an adiabatic MSW-transition between active and sterile neutrinos takes place, it 
could dramatically diminish the electron neutrino flux. 
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1 Introduction 



Sterile neutrinos closely degenerate with active neutrinos may have, if they exist, escaped 
detection in the laboratory and astrophysical experiments performed so far. They can, however, 
reveal themselves through measurable oscillation effects in astronomical-scale baselines P]Gj- 
They may affect the relative fluxes of the active neutrinos at the Earth if the mass-squared 
difference obey 5m 2 > E/L, where L is the distance to the source and E is the energy of 
neutrinos. With ultrahigh-energy cosmic ray (UHECR) neutrinos one can reach the sensitivity 
of Sm 2 > 1CT 18 eV 2 in the future neutrino telescopes like the ICECUBE PJ. It should be 
emphazised that the mass-squared differences below 5m 2 < 10 -11 eV 2 cannot be probed in 
solar neutrino, atmospheric neutrino or laboratory experiments. The most stringent present 
constraint on the active-sterile neutrino mixing comes from cosmology |<5m 2 | sin 2 2ip < 
5 x 10 -8 eV 2 . (These cosmological bounds may be avoided in the case of large lepton number 
asymmetries in the early universe |5J). It is generally thought that neutrinos are produced in 
UHECR sources via a pion-muon decay chain, which yields the flux ratios F® : F® : F® = 
1:2:0, where F® is the the flux of the neutrino flavour v a . If one considers only the active 
neutrinos v e , and v T and takes into account their "bi-large" mixing behaviour, observed in 
the solar and atmospheric neutrino measurements, the flavour ratios of neutrinos at the Earth 
are predicted to be F e : F^ : F T = 1 : 1 : 1 0. In [I] we found that active-sterile mixings 
with a tiny 5m 2 can change these ratios at the level of tens of percents. The effects of this size 
should be easily detected in the future experiments. 

In the present paper we shall extend our analysis of the degenerate active-sterile mixing 
to supernova neutrinos. As far as oscillations are concerned, the situation is for supernova 
neutrinos quite different from that for the UHECR neutrinos. In the case of UHECR neutrinos, 
the vacuum oscillations play a central role while for the supernova neutrinos matter effects are 
decisive and vacuum oscillations have usually no effects. In the case of supernova neutrinos the 
energy is also much lower, so they could be studied with experiments like HyperK and UNO. 

In the dense core of the supernova the neutrino Hamiltonian is extremely matter dominated. 
The interaction eigenstates, in which neutrinos are produced in various weak interaction pro- 
cesses in the core, coincide in a good accuracy with the eigenstates of this Hamiltonian. In 
transit through the envelope to the surface of the star, the flavour composition of these eigen- 
states changes. Neutrinos also pass through the MS W resonance regions corresponding to the 
solar and atmospheric oscillation scales, which may affect their behaviour. When leaving the 
star, neutrinos are in mass eigenstates that consist of different flavours according to the mixing 
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pattern they have in vacuum, and they will propagate as those states to the Earth without 
further oscillations. Consequently, if sterile neutrinos exist and if they mix with the active 
neutrinos, they will be present in the mass eigenstates entering the Earth. They will thereby 
affect the fluxes of the active neutrinos measured in neutrino detectors. 

There is an interesting possibility that one or more of the degenerate active-sterile neutrino 
pairs encounter a resonant mixing in the outer skirts of the envelope of the star. This is 
possible if that pair has \5vn?\ > 10~ n eV 2 , which corresponds to the oscillation length of the 
order of the giant progenitor star radius. Whether or not such a resonant mixing really occurs 
depends on how well the adiabaticity conditions are met in the resonance region, which in turn 
depends on the details of the density profile of the envelope in its outer skirts. If it does occur, 
it may change the flux ratios of neutrinos and antineutrinos dramatically, in particular if the 
active-sterile vacuum mixing angle is small. 



2 A sterile mixing scenario 

In the three- flavour framework the flavour fields v e ,v^, v T and the mass eigenfields vi,v~2, ^3 
are related to each other through v\ = Uui>i , where U is a mixing matrix, in the following 
parametrized as 

C12C13 S12C13 Sl3 ' 

U = -S12C23 - C12S23S13 C12C23 - S12S23S13 S23C13 j (!) 
\ S12S23 - C12C23S13 -C12S23 - S12C23S13 C23C13 J 

where Cjk = cos Ojk and Sjk = sin (9^ (we neglect the possible CP violation). The solar neutrino 
data, together with the recent KamLAND reactor data [Zj, indicate the mixing angle #12 to 
be bounded into the range 0.50 < 6*12 < 0.67 (Hj, On the other hand, observations of 
atmospheric neutrinos are consistent with maximal mixing between the mass eigenstates vi 
and U3, their mixing angle being within the range 0.64 < #23 < 0.96 [9], [10J. The third mixing 
angle, #13, is bounded by the results of CHOOZ and Palo Verde fTJJ to small values, 
< 013 < 0.1. 

Let us now assume that there exists three sterile neutrinos u s \, v S 2 and v S 3, which mix 
pairwise with the states i>i, z>2 and £3. We denote the new mass eigenstates that result from 
this mixing as follows (i = 1, 2, 3): 

Vi = cos Lfi Vi - sin w v si , 
v i = sin tpi Vi + cos ipi v si . 
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The antineutrino states are defined similarly. We assume that the mass difference of the states 
Ui and v[ is so small that in the processes, like particle decays, which are measured in laboratory 
experiments, these two states are not distinguished but appear as a single state with couplings 
equal to those of the active state V{. There are models where degenerate neutrino pairs may 
naturally appear, see e.g. ^Hj- It should be emphasized that in some of these models like in 
the well-known Pseudo-Dirac model, the mixing angles are naturally close to their maximal 
value of j. 

In the presence of the sterile neutrinos, the neutrino mixing matrix Q is modified to the 
form 



' COS (fi U e l COS (f2 U e 2 COS <£>3 U e s sill (pi U e \ sill (p% U e 2 sill 933 C/ e 3 ' 

cos (pi U^i cos tp2 U u2 cos (p 3 sin (p t U u i sin <p 2 U u2 sin (p 3 U ^3 

cos (pi U T i cos ip2 U T 2 cos <ps U T % sin ipi U T i sin if2 U T 2 sin 993 U T s 

— sin tpi cos ipi 

— sin (p2 cos (f2 

— sin ipz cos (^3 



Obviously, this is just one possible way to realizise the mixing between three active and three 
sterile neutrinos, not the most general case. 



3 Fluxes of supernova neutrinos 

We shall now study the effects of the sterile neutrinos on the fluxes of supernova neutrinos. 
As mentioned, for supernova neutrinos matter effects play the main role, in contrast with the 
case of UHECR neutrinos, where vacuum oscillations are important. Because of the effects of 
matter, the flavour composition of the supernova neutrino flux observed at the Earth differs 
from that in the production region (see, e.g. [T^ ITH]). 

The neutrinos and antineutrinos may undergo matter enhanced MSW transitions inside 
the star if 5m? < 10 4 eV 2 . In the case of three active neutrinos there are two possible MSW- 
resonance regions, at the densities 

pu ~ 10 3 - 10 4 g/cm 3 , p L ~ 10 - 30 g/cm 3 . (4) 

The subscript H refers to the so-called high-resonance region, which corresponds to the atmo- 
spheric neutrino oscillation (<5m 2 tm = 2.6 x 10~ 3 eV 2 ,sin 2 # atm = 0.52 ^9j), and the subscript L 
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refers to the low-resonance region, which corresponds to the large mixing angle (LMA) solar 
neutrino oscillations (Sm^ = 6.9 x eV 2 and tan 2 6 Q = 0.43 !.)]). In the case of the normal 
mass hierarchy (mi < 777,2 <C 7713) both resonances occur for neutrinos, whereas in the case 
of the inverse mass hierarchy (7723 <C m\ < 777-2) the high-resonance occurs for antineutrinos 
and low-resonance for neutrinos. If the system is not fully adiabatic, the MSW effect is not 
complete but a level crossing from one matter eigenstate to another will occur in the resonance 
region. The level crossing probability, the so-called Landau-Zener probability, is given in [16J. 

The matter effects depend on the density profile of the progenitor star. It can be shown 
(e.g. [151) that 

propagation through the low-resonance region, determined by the solar neutrino 
parameters, is adiabatic, i.e. the Landau-Zener probability for it is Pl = 0. For the high- 
resonance one can distinguish three cases, defined by the values of sin 2 #13 and neutrino energy 
and differing in the values of the Landau-Zener probability (Ph) 

1. Adiabaticity breaking region: sin 2 #13 < 1(T 6 x (E/10 MeV) 2 / 3 , where P H « 1; 

2. Transition region: sin 2 13 ~ (1(T 6 - 10 -4 ) x (E/10 MeV) 2 / 3 , where < P H < 1; 

3. Adiabatic region: sin 2 #13 > 1(T 4 x (E/10 MeV) 2 / 3 , where P H » 0. 

In the adiabaticity breaking region the H-resonance has no effect to the evolution of the neutrino 
states, whereas in the adiabatic region a full conversion will occur. For simplicity, we will 
concentrate in what follows on these two extreme cases 1 and 3, omitting the transition region 
case 2. 

In the standard case of three active neutrinos with the normal mass hierarchy the fluxes 
of the mass eigenstates leaving the star and entering the Earth later on, which we denote 
by F = (Fi, F2, F3), are obtained from the production fluxes of the flavour eigenstate in the 
supernova core, denoted by F° = (P e °, P° P T °), through jTH] 



describes the conversion probability inside the star. The fluxes of different flavour eigenstate 
neutrinos v a at the Earth are then given by 

3 



F = PF 1 



(5) 



where the matrix 



/ PhPl 
P= P h (1-Pl) 



1 - Pi (1 - P h )Pl ^ 
P L (l-P H )(l-P L ) 

Ph j 



(6) 




(7) 
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Let us now consider the situation in the case of three additional sterile neutrinos. Sterile 
neutrinos are not produced in the weak interaction processes taking place in the core of a 
supernova. Therefore, the fluxes of the mass eigenstates on the surface of the supernova are, 
like in the standard case, determined by the production rates of the active neutrinos in the 
core. Nevertheless, during the propagation of neutrinos through the envelope of the supernova, 
sterile components ^ s i)^s2,^s3 are developed as a result of the active-sterile mixing, that is, 
active-to-sterile transitions take place. The flavour composition of the neutrino flux on the 
surface and at the Earth is given by 

6 3 3 

F <* =J2\ U ^\ 2Fi =I>osVilM 2 ^ + ^sinVl^'| 2 *V (8) 

t=l i=l j'=l 

where a = e,fjt, r, si, s2, s3 and F^s are given in Eq. JH} instead that P is now six-dimensional 
and F° = F°, F T °, 0, 0, 0). If the active-sterile mixing angles cpi are all equal, the relative 
fluxes of different flavours, F a , do not differ from those of the standard case, F^ M . In partic- 
ular, this is true in the case where all the mixing angles (fi are close to their maximal value of 
?, as predicted by some models. 

In the case of the normal mass hierarchy, antineutrinos do not encounter MSW resonant 
conversions, and they end up to different mass eigenstates than the corresponding neutrinos. 
The counterpart of the matrix P for antineutrinos is a unit matrix. Hence the active-sterile 
neutrino mixing affects antineutrinos and neutrinos differently, and consequently the ratio 
Fa/ Fa generally differs from its value in the nonsterile case. 

In the case of the inverted mass hierarchy, the high-resonance is encountered by antineutri- 
nos and the low-resonance by neutrinos. Obviously, the fluxes of neutrino mass eigenstates on 
the surface of the star are obtained from Eq. © by setting Ph = 1. It is also straightforward 
to see that the antineutrino counterpart of the matrix P given in Eq. JOJ, denoted by P, is 
obtained in the inverted mass hierarchy case from the matrix P by replacing Ph with Ph 
(actually Ph = Ph [Hj) and Pl with 1 — Pl ^H]- The fluxes of the active antineutrinos at 
the Earth are then given by 

3 3 

F a = ^cos 2 <pi\Uai\ 2 Fi + ^ sin 2 <Pi'\U ai > | 2 (9) 
i=i i'=i 

where -F = (F®, F9, F®, 0, 0, 0) are the fluxes of antineutrino flavours in the production region 
and F/s on the right-hand side are obtained by 

F = PF°. (10) 



6 



Although the low-resonance point itself does not appear in the antineutrino sector, a conversion 
effect is possible at the low-density region, and Pl takes into account the adiabaticity of this 
transition In our calculations we will assume Pl=0. 

4 Results 

In the following we shall present numerical estimates for the flux ratios of supernova neutrinos 
at the Earth in the case of degenerate sterile neutrinos, and compare them with the corre- 
sponding results in the standard, nonsterile case. We allow the ordinary mixing angles 0^ vary 
in their phenomenologically allowed regions quoted in Section 2. The active-sterile mixing 
angles cpi are allowed to vary arbitrarily in the range from ifi = (no mixing) to ipi = 7r/4 
(maximal mixing). 

There exists a considerable uncertainty concerning the initial fluxes of neutrino flavours 
in the production region in the core of the supernova. The value of the ratio F e ° : F g ° : 
where F® = F® = = F® = F®, varies according to the model one uses Traditionally 
these fluxes are supposed to be equal, i.e. F° : F® : F® = 1 : 1 : 1, which is one reference 
value we will use in our analysis. According to recent detailed studies of microprocesses taking 
place in supernova core, the flux ratios may considerably differ from this simple assumption. 
From the results of [20J we obtain, after integrating over the energy spectrum, the ratios 
F® : Fg : F® = 4 : 3 : 2, which we will use in the following. 

Let us first assume that no active-sterile matter conversion takes place in the outer skirts of 
the star. The opposite case will be considered in the end of this chapter. We compute the flux 
ratios F e /F a , F e jF^, and F g /F a , where F a is the sum of the fluxes of active neutrinos other 
than v e and v e . The results are shown in tables^andEl In table^we have compared the values 
of relative fluxes with and without active-sterile vacuum mixing, for the normal and inverted 
mass hierarchy and for an adiabatic and non-adiabatic conversion at the high-resonance. We 
see immediately, that in the case of non-adiabatic transition in the high-resonance the ratios 
are independent of the mass hierarchy. This is in agreement with the conclusion of ^E], that 
for sin 2 #13 < 10~ 6 the observation of supernova neutrinos are insensitive to the mass hierarchy. 
In adiabatic case, the relative amount of electon neutrinos is larger for the inverted hierarchy. 
The mixing with sterile neutrinos widens the range of the possible values of the flux ratios in 
all cases. With equal initial fluxes the situation is more simple, as shown in table in the 
sense that the ratios are independent of the neutrino mass hierarchy and the adiabaticity of 
the MSW conversions. Also in this case the active-sterile mixings considerably widen the range 
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of the possible flux ratios. 







Normal hierarchy 


Inverted hierarchy 






active 


active + sterile 


active 


active + sterile 




F e 
F a 


0.18 


0.12 - 0.28 


0.22-0.25 


0.16-0.33 


adiab. 


Fe 
Fa 


0.28 - 0.29 


0.18 - 0.46 


0.20 


0.14 - 0.31 




F e 
Fe 


0.63 - 0.67 


0.61 - 0.72 


1.07-1.21 


0.99 - 1.35 




Fe 
Fa 


0.24-0.27 


0.17-0.33 


0.24-0.27 


0.17-0.33 


non-adiab. 


Es. 

Fa 


0.30 


0.20 - 0.47 


0.30 


0.20 - 0.47 




Fe 
Fe 


0.77 - 0.92 


0.69 - 1.10 


0.77 - 0.92 


0.69 - 1.10 



Table 1: Results with initial flux ratios 4:3:2 



In a summary, the effect of the active-sterile mixing is in general less striking for the 
supernova neutrinos than what we found in our previous study flj it to be for the UHECR 
neutrinos. The relative flux of the electron antineutrino can, however, differ from its nonsterile 
value as much as 50 % if the active-sterile mixing angles (fi for different degenerate pairs 
differ suitably from each other. The basic reason for this is that the electron neutrino and 
antineutrino end up to different mass states in the surface of the star, so that their fluxes are 
sensitive to different active-sterile mixing angles. 





active 


active + sterile 


Fe 
F a 


0.25 


0.17-0.38 


Fe 
F a 


0.25 


0.17-0.38 


Fe 
Fe 


1.00 


1.00 



Table 2: Results with initial flux ratios 1:1:1. These values are valid for both normal and 
inverted mass hierarchy and independently of adiabaticity. 

So far we have assumed that the MS W conversions between the active and sterile states 
do not occur, in other words, the transitions between the mass states V{ = cosi^j V{ — sin<^j v s i 
and v[ = sin (fi Oi + cos ifi v s i are non-adiabatic. The adiabaticity depends on the profile of 
the envelope at the resonance region, as well as on the mixing angles cpi and the mass-squared 
difference mf — mf. In order to get a feeling of what may happen, we have studied, utilizing 
the analysis of ref. |21j . a case of mixing of one degenerate active-sterile neutrino pair more 
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closely, and we have found that for mf — mf < 10~ n eV 2 the transition is highly non-adiabatic 
for mixing angles up to about 35° and adiabatic only in the case of nearly maximal mixing. 
For larger mass-squared differencies adiabaticity can be reached with smaller mixing angles. In 
general there will be a range in the parameter values, where the transiton is partially adiabatic. 
Obviously, our analysis does not make full justice to the problem, but it does convince us that 
the sterile mixing of the sort we have been looking at may have detectable effects. 

Let us now assume that these transitions are fully adiabatic and the transitions v; L — > v[ 
do take place. In the case the mixing angle ipi is close to maximal this would not have 
any influence on the fluxes of the active neutrinos. In contrast, if the mixing angle (fi is 
small, the transition would mean a conversion of a predominantly active neutrino state into a 
predominantly sterile state, which would diminish the observable flux of supernova neutrinos 
at the Earth. The mixing angle (pi cannot be, however, arbitrarily small, since for very small 
angles the transition would not be adiabatic (unless the mass-squared difference is large). Large 
effects are anyway possible, and if the mixing angles of different degenerate pairs differ from 
each other, the flux ratios may change dramatically. This can be seen in Table 3, where the 
flux ratios are presented for the case of the normal mass hierarchy and assuming adiabatic 
Vi — > v[ transitions. The mixing angle ifi is allowed to vary in the range 5° — 45°. The relative 
flux of the electron neutrino can substatially decrease as a result of these conversions. 







Normal hierarchy 






F e ° : F g ° : F° = 4 : 3 : 2 


F° : Fl : F* = 1 : 1 : 1 




F e 
F a 


0.004 - 0.21 


0.005 - 0.24 


adiab. 


Fe 
Fa 


0.30 - 0.64 


0.26 - 0.49 




F e 
Fe 


0.006 - 0.64 


0.009 - 0.97 




Fe 
F a 


0.005 - 0.25 


0.004 - 0.24 


non-adiab. 


Fe 
F a 


0.30 - 0.66 


0.26 - 0.49 




Fe 

F e 


0.009 - 0.88 


0.008 - 0.97 



Table 3: Results with adiabatic fi — > v[ -transitions 



In considering the significance of the effects discussed above one should also consider the 
anticipated accuracy of the experimental determination of the supernova neutrino fluxes. For 
example, in the SNO detector, for a typical supernova in the Galaxy with a distance of 10 kpc, 
the relative uncertainty of the v e and v a fluxies is about 5 % and that of the v e flux about 10 % 
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|22j . Consequently, the uncertainties of the determination of the flux ratios are roughly 5 % for 
F s /F a and roughly 10 %-15 % for F e /F a and F e /F s . In the Super-Kamiokande detector the 
i> e flux can be determined more accurately, with the uncertainty of about 1 % (corresponding 
to about observed 10 000 events). The accuracy of the v a flux determination is estimated to 
be about 4 % (about 700 events) (221- The v e flux is, in contrast, quite hard to determine at 
the Super-K. While noting that the effects of the active-sterile mixing we have discussed would 
be detectable in these detectors, we must stress the importance of developing methods of an 
accurate determination of the flux of electron neutrinos in future neutrino experiments. 

5 Summary 

We have investigated how supernova neutrino fluxes are affected by the existence of sterile 
neutrinos closely degenerate with active neutrinos. Sterile neutrinos are not produced in the 
core of a supernova, but in the matter eigenstates entering the surface of the supernova envelope 
sterile components are developed. We have used two sets of initial fluxes of neutrinos in the 
production region, F e ° : F g ° : F° = 1 : 1 : 1, and F e ° : F g : F° = 4 : 3 : 2. Matter effects 
were taken into account and we have explored both normal and inverted mass hierarchies. 
The effects caused by the degenerate sterile neutrinos will be detectable in future neutrino 
experiments. Particularly large effects are possible for the electron antineutrino flux. 

It is possible, at least in principle, that the degenerate mass eigenstate pairs encounter a 
MSW resonance conversion in the outer skirts of the supernova envelope. This may dramat- 
ically change the relative fluxes of the neutrinos interacting in detectors. If the active-sterile 
mixing angles is small, a transition from a predominantly active state to a predominantly 
sterile state will occur, which can diminish e.g. the electron neutrino flux a lot. 
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